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Abstract. For analyzing anisotropic low relative- velocity correlation- functions 
^ ' and the associated emission sources, we propose an expansion in terms of carte- 

ls , sian spherical harmonics. The expansion coefficients represent angular mo- 

' ments of the investigated functions. The respective coefficients for the correla- 

tion and source are directly related to each other via one-dimensional integral 
transforms. The shape features of the source may be partly read off from the 
ly^ , respective features of the correlation function and can be, otherwise, imaged. 
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1. Introduction 



. , Use of interferometry for the determination of reaction geometry came to nuclear 

■ physics from astronomy. In both fields, the intensity technique has been first used 

for determining gross sizes of the radiating regions. While astronomy has since 
moved on to a systematic model-independent determination of the fine details of 
observed objects within the phase version of the method, the nuclear physics has 
not advanced beyond fitting the experimental data with models that parameterize 
coarse shape- features of particle-emitting regions [E| O]. Imprinted in the shape 
features can be the reaction history and, specifically, differences in emission times 
for different particles and stalling of a system when passing through the quark- 
gluon phase transition. This talk is dedicated to the systematic determination of 
the shape features of emitting regions in reactions with multi-particle final states. 
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2. Low- Velocity Correlations in Reactions 

Ability to learn on geometry of emitting regions in reactions relies on the possibility 
of factorizing the amplitude for the reaction into a wavefunction <i>q^'' for the pair 
of detected particles and an amplitude remnant. The wavefunction can be com- 
puted and, at at low relative velocity, its square may exhibit pronounced spatial 
features such as associated with identity interference, resonances or Coulomb repul- 
sion. Those features can be regulated by changing the relative particle momentum q. 
The amplitude remnant, squared within a cross section and summed over the un- 
observed particles and integrated over their momenta, yields S' which represents 
general features of reaction geometry, without a significant variation with q. When 
looking at the inclusive two-particle cross-section then, the structures in |$(+) p can 
be used to explore the geometry in S': 



The naive expectation, met at large q for a multi-particle final state, is that the 
emission of two particles is uncorrelated. Thus, one normalizes the two-particle cross 
section with a product of single-particle cross sections, and looks for the deviations 
of the thus constructed correlation function C from 1 : 



7^(q) = C(q) - 1 = {^'["l^ - 1 = / dr (|$(-)(r)p - l) 5p(r) . (2) 



The last equality follows from the fact that, at large g, the relative wavefunction 
squared is equal, on the average, to 1; in combining this with C ~ 1 at large q, the 
source S, following the cross-section normalization, turns out to be normalized to 1, 
/ dr5'(r) = 1. Equation Q links the deviations of correlation function C from 1, 
to the interplay of deviations of the relative wavefunction from 1 with the geometry 
in S. With S normalized to 1, it may be interpreted as the probability of emitting 
the two particles at the relative separation r within the particle center of mass. The 
emission is integrated over time, as far as S is concerned. 

3. Analysis of Correlations 

Knowing 7?.(q) and |<i>|, on can try to learn about S; mathematically, this repre- 
sents a difficult problem involving an inversion of the integral kernel i^(q, r) = 
|$q~''(r)p — 1, in three dimensions. The situation is simplified by the fact that K 
depends only on q, r and the angle ^qr, and, correspondingly, can be expanded 
in Legendre polynomials, i^(q, r) = + I) Ke{q,r) {cos 9). If the corre- 

lation and source functions are expanded in spherical harmonics Y^™, Ti{q) = 
^/Att X^fm T^^"^ (<?) Y^™ (q) : one finds that the three-dimensional relation Q is equiv- 
alent to a set of one-dimensional relations for the harmonic coefficients [ 4 




(1) 



1 



da 




(3) 



Angular Moment Analysis of Correlations 



3 



For weak anisotropics, only \ow-£ coefficients of 7?. or 5 are expected to be significant. 
The £ = version of Q connects the angle-averaged functions. 

The above suggests a systematic analysis of the correlation functions and source 
functions in terms of the harmonic coefficients. A problem, however, is that it is 
cumbersome to analyze real functions in terms of complex coefficients lacking a clear 
interpretation for £ > 1. This suggests looking for another basis for the directional 
decomposition of correlation functions and sources, that would be real and have a 
a clear geometric meaning. Such a basis may be constructed starting from a unit 
direction vector fia = (sin 9 cos 0, sin 6 sin 0, cos 9) . 

The tensor product of £ vectors fia yields a symmetric rank-£ cartesian tensor 
that is a combination of spherical tensors of rank £' < £ and the same evenness as £: 

)ai...ai = ^ai ■ ■ ■ — ^ ^ ^£'ni ^ • (4) 

£'<e,m 

A projection operator 2? may be constructed in the space of symmetric rank-^ 
tensors, out of a combination of Kronecker (5-symbols, that makes a symmetric 
tensor traceless, 

Xl(2?nOaaa3...a, =0. (5) 

a 

The tracelessness of iVn^) ensures that the products (Pn^)^^.. ^^ are solutions 
of the Laplace equation and, thus, the {T>n^)a-^,,,ai are combinations of spherical 
harmonics of rank £ only. The components are real and may be used to replace Y^"^ . 
The lowest-rank tensors are: 

{T)fl )q,^ Q3 — "^"(^qi a2 ^'^03 ~^ ^ai Q3 ^a2 ~^ ^ck2 ct3 ^cti) ■ (6) 

5 

The completeness relation in terms of the cartesian components is 

I 0.1. ..a.1 

1 (2£-f 1)!! 

= J\ 2^ "ai ••■'^a, , (7) 

^ OL\...ai 

where the second equality follows from V = T)^ = . The completeness relation 
can be used for expanding 7?, or 5* in terms of cartesian tensor components 

7^(q)= Un'5{n'-^)n{c()^Y. E (8) 

where the coefficients are angular moments, 7^ai...Q^ (g) — / 7^(q) (X>g^)ai...a^. 
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With the TZ and S cartesian coefficients being identical combinations of the 
respective y^™-coefficients, the corresponding TZ and S cartesian coefficients are 
directly related to each other, 

TeL^^.a, (g) = 4^ y dr K,{q, r) s'^.-o,, (r) . (9) 

For weak anisotropics, only low ^-coefficients for the source or correlation matter, 

7^(q) = 7^(") {q) + 7^(:) (g) go + E ^"lU (?) 9-2 + ■ • ■ • (10) 

Here, 7?.(°^ is the angle-averaged function. The dipole and quadrupole distortions, 
T?.^^^ and TZ^'^\ can be represented in terms of amplitudes and distortion vectors: 

7^« = e« and 7^S(,) = ^ efj eg e^j eg - (i?g + i?g) e^^^ eg, 

respectively. 



4. Illustration 



We illustrate the analysis in terms of cartesian harmonics taking an anisotropic 
Gaussian source, elongated along the beam axis, and displaced along the total pair 
momentum P, at an angle of 30° relative to the beam axis, cf. Fig. ^ Figure |21 

Displaced deformed source 

S(x,0,z) [fm-'i 




5 beam ' 



moments evaluated moving 
at constant r about the origin 

Fig. 1. Sample Gaussian emission source within the plane of relative directions, 
with the z-axis pointing along the beam and the a;-axis along P^. 



shows low-^ source values and directions, from the sample source decomposition in 
terms of the cartesian harmonics. The angle-average source maximizes at r = 0, 
although the original source is displaced from the origin. The distortions can only 
become significant at distances comparable to source size; in fact, at small distances 
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Fig. 2. Low-^ values (left) and angles (right) for the relative source, as a function of 
the relative distance r. Lines and symbols represent, respectively, the exact values 
and values obtained through imaging of the anisotropic correlation function. 



For the anisotropic source, we next construct a correlation function, assuming 
the classical limit of repulsive Coulomb interactions, suitable for emitted interme- 
diate mass fragments. In that limit, the kernel K \s a. function of of 0qr and t/tc, 
where rc is the distance of closest approach in a head-on collision, ^ — — ^° ^ , 
as 

l^|2 _ d^go _ e (1 + cos6lqr - 2rrjr) (1 + cos6>qr - r^/r) ^^^^ 
"^^^ y^(l + cos V)^ - (1 + cos V) 2rc/r 



The 1 = kernel is Kq — 0(r — rc) ^1 — rc/r — 1, while the higher-rank kernels 
may be obtained numerically. 

In the classical limit, the correlation function reflects the distribution of relative 
trajectories emerging from an anisotropic source, cf. Fig.|3| In the low-energy limit 
of r — Tc ^ Tc, the trajectories turn back away from r = 0, lumping around the 
direction of r they emerge from. For r,. comparable to the source size, the low- 
energy limit holds and only the source margins contribute to the correlation. In 
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the high-energy limit of ^ r, the trajectories emerge isotropically, except for the 
shadow left by the r < region in the — r direction. For small compared to the 
source size, the high-energy limit generally holds and the correlation represents the 
integral features of the source, as nearly all source points contribute. 




Fig. 3. A schematic relative source with repulsive Coulomb trajectories super- 
imposed coming out isotropically from position r within the source. The region of 
radius Vc is inaccessible to the trajectories. 

Characteristics of the classical Coulomb correlation function for the source in 
Figs. Hand 121 are next shown in Fig. 01 plotted there vs. r^^^^ oc q. The dipole angle 
for C follows the variation of the dipole angle for S. On the other, the quadrupole 

— 1/2 

angle exhibits a jump by 90° as a function of rc . This jump is associated with 
the change in sign for K2 and the resulting prolate-to-oblate shape transition for C . 
For more schematic sources, one or more low-£ correlation amplitudes may vanish 
and/or angles may exhibit less variation. 

We further test out the imaging [|2 IHI of the features of our anisotropic source 
from the correlation function. We assume that the £ < 2 cartesian coefhcients have 

— 1/2 

been measured for our Coulomb correlation function, at 80 values of rc . We 
generate the coefficient values randomly assuming a constant r.m.s. error of 0.015 [ 

From the simulated measurements, we restore the coefficients for the source, 
within the region of r < 17 fm, at the number of points, in r, that drops with t. 
5 for f = 0, 4 for I = \ and 3 for £ = 2. The results of restoration, represented by 
symbols in Fig. [3 compare reasonably with the original source features. 

Given the cartesian coefficients for the source as a function of r, it is straightfor- 
ward to compute the cartesian moments of the source. The moments from imaging 
are compared in Table^to those computed directly and they are found to reproduce 
the latter reasonably. 
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0.00 0.25 0.50 0.75 1.00 1.25 

r„"'^^ (fin"'^^)'«-inverse square root 
of Coulomb radius 

Fig. 4. Low-£ characteristics of the Coulomb correlation function for our sample 
source: amplitudes (top panel) and angles (bottom panel) as a function of the 
inverse square root of the head-on return radius Tc- 



5. Conclusions 

Characteristics of the anisotropic correlation functions in reactions, and of corre- 
sponding emission sources, may be quantitatively expressed in terms of the cartesian 
coefficients that represent angular moments of the functions. The coefhcients for 
the correlation and the source depend on magnitude of the relative momentum and 
separation, respectively. The respective coefficients for the source and correlation 
are related to each other through one-dimensional integral transforms. Certain 
shape features for the source may be directly read off from the correlation features 
in terms of the coefficients. Otherwise, the source features may be imaged. 
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Table 1. Cartesian moments and weight of the source at r < 17 fm, obtained in 



imaging and by direct integration. 

tjuantity 



47r/drr2 5W 

(x) 

((x-(x))2)l/2 

{{z-{z))Y" 
{x)){z-{z))) 



Unit 



fm 
fm 
fm 
fm 
fm 
fm^ 



Restored Original 



0.99±0.05 
2.47±0.11 
4.25±0.13 
3.80±0.24 
3.81±0.22 
5.54±0.19 
2.23±1.49 




